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^ ■ Abstract 

We demonstrate that an effect other than anharmonicity can severely distort the spectroscopic 
signatures of quantum mechanical systems. This is done through an analytic calculation of the 
spectroscopic response of a simple system, a charged torsional pendulum. One may look for these 
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effects in the optical data of real systems when for example a significant rocking component of rigid 

c/j , polyhedra plays a significant role in the lattice dynamics. 
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In a typical spectroscopic experiment, electromagnetic radiation couples to the charge 
degrees of freedom of a system whose underlying dynamical behavior is described by a 
Hamiltonian. In many cases, the behavior of a system near a stable equilibrium can be 
adequately described by a simple harmonic Hamiltonian, a choice that is often well justified 
as an expansion in some suitably chosen coordinate about a minimum of a more compli- 
cated potential. The potential expansion coordinates are often chosen to be one or more 
of the Cartesian coordinates x, y, and z, such as is done in the quantum theory of lattice 
vibrations 

The coupling of radiation to matter also involves matrix elements of the Cartesian vari- 
ables and as a result, selection rules arise which forbid optical transitions between vibrational 
levels which are not adjacent in energy. The effect on spectroscopy is to produce a single 
peak in the system response at the oscillator's fundamental frequency. The symmetry for- 
bidding the transitions to higher energy levels is imposed by the initial choice of Cartesian 
expansion coordinate, a choice that in certain situations could be improved upon toward de- 
scribing the system dynamics. An example of the latter is a system with inherent curvilinear 
geometrical constraints. 

In this manuscript, we show that for a simple and familiar model system where curvilinear 
motion is inherent, profound effects on the spectroscopic response functions are realized. We 
exemplify this principle by studying the dynamical response of a torsional pendulum with 
charge degrees of freedom. 

The quantum mechanical problem which is analogous to the classical torsional pendulum 
is shown in Figured). This situation corresponds to a situation where 360° rotations of the 
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particle wavefunction produces a distinct state, which cannot interfere with the unrotated 
state. The "particle on a pig's tail" system shown in Figure ^ is can be viewed as a ID 
harmonic oscillator wrapped many times around a cylinder. We proceed to investigate the 
spectroscopic response functions of this harmonic system. 

***Quantum Lorentz Oscillator*** To begin we show that, in the absence of damping, a 
charged linear harmonic oscillator gives rise to response only at the oscillator resonant fre- 
quency and hence will have a single peak in the spectrum at this frequency. The Hamiltonian 
for the linear oscillator is 

rr P* 1 2 2 

where m is the particle mass and Uq is the resonant frequency of the pendulum. The 
momentum operator p is conjugate to the position coordinate x, which can have any value 
on the x axis. We exploit extensively Dirac's factorization procedure in order to find the 
current-current correlation function, and hence the linear response of this model system. 
The Hamiltonian in quantized form is: 

H = (a) a + -)Hujo 

and 

H\n) = hu (n + ^)\n) 

Where \n) is the nth harmonic oscillator state and a) and a are the raising and lowering 
operators of Dirac's theory, which are related to the observables x and p by 



•r = \ / - - (a + a ] ) 

Zmujn 



mLd h t 

The spectrum of this system is a ladder of states in energy separated by Hwq. 
The vector components of the polarization operator are 

Px QX, Py 

and the associated current operator can be determined from the Heisenberg equation of 
motion, 

jx = ^[H,P X } = %. 
n m 
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The matrix elements we will be interested in are 



WAN = £(-i^)<»l(« -«*)!*.) 



Q . mu h 

= ){Pn,m-l - Vm+lJ 

m V 2 

The optical conductivity, cri(c<j), which represents the amount of current that is generated 
by an oscillatory field of frequency u, is generally related to these matrix elements through 

= ^- y2\(n\i\0)\ 2 (S(^-E n + E ) + 5(hu + E n -E )), 
Vuj / — J 

which for the case of the linear harmonic oscillator, is 

7T Q 2 

= V2^ - ^o) + + w )). 

The optical conductivity has a single peak corresponding to absorption of one quantum of 
radiation of frequency uq. When damping is included, this expression generalizes to the 
Lorentzian line shape, which is encountered in the context of the power absorption in the 
damped driven harmonic oscillator (2| and is commonly used in the analysis of optical data 
where a sharp resonance characterizes the absorption. 

The absorption of radiation at only a single frequency is a bit surprising since the har- 
monic oscillator has energy levels for n = 0, 1, 2, 3,..., and so one may expect that there be 
corresponding optical transitions with frequencies uj = (E n — E )/h = u , 2u , 3co , .... We 
shall see below that these transitions are forbidden for the linear harmonic oscillator but 
need not be in general. 

We now repeat this treatment for the torsional oscillator, shown in Figure ^p. One can 
view this as a ID harmonic oscillator wrapped on the surface of a cylinder, with the position 
variable x closely related to the angular coordinate 9. The Hamiltonian for the torsional 
oscillator is 



i? i 
— h -- 

21 2 

where / is the moment of inertia and ujq is the resonant frequency of the pendulum. The 
angular momentum operator L is conjugate to the (unbounded) angular coordinate 9, and 
the Hamiltonian describes a harmonic oscillator. We can again exploit extensively Dirac's 
factorization procedure as before with the appropriate substitutions x — > 9 and p — >• L. 
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Repeating the procedure outlined above for the linear harmonic oscillator, the vector 
components of the polarization operator are 

Px — Q r cos # 

P y = Qr sin 9 

and the associated current operators can be determined from the Heisenberg equation of 
motion, 

j x = ^[H,P x ] = -^(Lsme + sm9L) 
3y = ^[H,P y ] = ^(Lcos9 + cos9L). 

In the torsional oscillator, we are concerned with taking matrix elements of a nontrivial 
function of the angle 9, which appears quadratically in the Hamiltonian. This nonlinear 
dependence of the matrix element on the harmonic degree of freedom will permit transitions 
between states not adjacent in energy, an effect which is studied widely and has been termed 
electrical anharmonicity. The main effect of this for spectroscopy is that the conductivity 
will be a series of peaks at integral multiples of the resonant frequency We will demonstrate 
this effect below by calculating the spectroscopic matrix elements for this system. 

The matrix elements we will be interested in are of the form 

Qr 

(n\j x \m) = — — (n\(Lsm9 + sm9L)\m). 
Cribbing the result from the harmonic oscillator theory, 



iujoh, t . 

<a — a ] 
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we can write the current matrix elements in terms of matrix elements of the trigonometric 
functions of 9: 



(n\j x \m) = — ^—)(n\((a — a 1 ) siia.9 + sin9(a — a^))\ni) 



i-^y\ ~^~ (( n + l\Vn + 1 sin#|m) — (n — llv^sin^lm) 

_ ./ V — 

+{n\ sin 9\pm\m — 1)) — (n\ sin 9^/m +T|m + 1) 
.Qr y ^^(§ n+l mV / n + i _ E> n ^ hm y/n 

+§ ni „ t _iv / m - E> n<m+ Wm+ 1) 



= l 2I 



where 

§n,m = { n \ sm9\m) 
A nearly identical calculation for the current component j y gives 

(n\j y \m) = -i^J~^^(C n+1>m Vn + 1 - C n _i i , nV / ^ 

+C„ im _i V / ^ - C n , m+1 Vm + 1) 

C n ,m = (w| cos#|m) 

To take a general approach, we will calculate the matrix elements of the operator 

E n ,m = (n\e id \m) 

and use them to determine those of the trigonometric functions through deMoivre's identity. 

Again using a result from the harmonic oscillator theory, and introducing the factor c for 
notational simplification, 





E n<m (c) = (n\e^ a+a ^\m). 



The trigonometric matrix elements are then 

Sn,m(c) = ^(E n>m (c) - E rl)TO (-c)) 

C„, m (c) = ^(E n>m (c) + E n , m (-c)) 

To simplify the calculation of E„ m (c), one can invoke the Baker- Campbell-Haussdorf 
theorem [2|: 

e A+B = e A e B e -[A,B]/2 

which holds provided that A and B both commute with their mutual commutator. This is 
true for a and a\ so 

E n , m = (n\e ica e icai e c2[a ' a ^ /2 \m) = (n\e ica e icat e c2/2 \m). 

We can now expand the exponentials and operate on the bra and ket multiple times with 
the a and a'. Using the identities 

\m) = 4=10) 
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and 

(n| = (0|^4 
vn 



(zcr 



E„, m = e c2 /2 (n | e -^^L( G t )fc | m) 

fc=0 



e c2 / 2 V(n|e ica ^(at) fc ^£|0) 



cV 2 V/ n ie-M v3^±M (£T^ |0) 

^ fc! v 7 ^! v/(m + fc)! 



e 

fc=0 



r 2 /2 / i i™ (^) fc \/( m + i , \ 

= e c/i > ne tca ^— - vv ; - \m + k) 

Similarly, the e* ca can be expanded to act on the bra (n\, 



OO 

„2 i 



v ml 



OO OO 



, - ■■■■> sr^ sr^ (ic) 1 , sl (ic) k J (m + fc)! . 

e 1 E Ex (n|(fl) fc! /^t |m + } 



fc=0 1=0 " v m ' 

/! i/til fc! V 777,1 



£Z£Z /! 1 V(» + 0« ^! fc! vm 
Combining factors and using the orthonormality of the harmonic oscillator states, 

» (zc) fc+ V(/ + n)!(fc + m)! ^ 
Vn!m! ^ /!fc! 

e c2 / 2 ^ (ic) 2k+m - n (k + m)\ 



Vn!m! ~ fc!(fc + m-n)! 

e c2 / 2 (zc) m " n y> (-c 2 ) fc (fc + m)\ 
\frim~\ k\{k + m - n)\ 

e c2/2 (ic) m - n 



S(c 2 , m, n) 



\Jn\m\ 

The sum S(c 2 ,m,n) is a real function and can be expressed in terms of hypergeometric 
and gamma functions^. This expression can now be used to give the trigonometric matrix 
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elements of 9 and these are 



and 



S n , m (c) = -==S(c 2 ,m,n)^ ^ 

Vn!m! ^ 

1 e c2//2 2 [0 m — n even 

= ~ /=j=? £( c 

i Vn!m! [ (ic) m " n m - n odd 



C n , m (c) = -==,S(c 2 ,m,n)^ ^ 

vn!m! 1 

e c2 / 2 2 f (^c) m ~ n m — n even 

^(c , m, n) • 



Vn\m\ ' [ o m - n odd 

We restrict our attention at this point to the case of zero temperature. In that case, we 
are interested in transitions from the ground state, so that n = 0: 



(0\j x \m) = i^jJ — ^(§i,m + So^-iv 7 ™ - S , m +i V«i + 1) 
This is zero when m is odd. For m even, 



V 17 ' 7 2/ V 2 ^ k\i l(jfe + m-l)!v^! 

(k + m- ly.^cY^^/m (k + m + l)\(ic) m+l ^/ m + 1 



(fc + m- l)!^(m- 1)! (fc + m + 1)! ^m+lj! 
y (~c 2 ) fc f 



.Qr c 2 /2 v^ (~c 2 ) fc f (fc + m)^)™- 1 (ic)™" 1 ™ (ic) m+1 

*2/V^r e k\i i T^t + v^y T^T 



.Qr /ic^fi c 2 /2 A (-c 2 ) fc (ic) m r (fc + 2m) c 2 • 



= i^M— e c2 / 2 



2/ V z — /,■:/ v 



{-(f)" {icy- ({K + im) C* "I 



.Qr (-<?)* (ic)™' 1 



We can now sum the series, which gives exponentials: 

k=0 

= lf\^ eC2/2 ^{- cVC2+e " C2 ( 2m + c2 )} 



'mi 

2/ V 2 V^T 
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FIG. 2: The matrix elements of the torsional oscillator for three values of angular uncertainty. Only 
the integral values of m are meaningful. Inset: The ground state wavefunctions corresponding to 
these three uncertainty values. 



The analysis of the (0\j y \m) matrix element follows similarly. The intensities of the 
transitions are determined from 

12 ml m odd 



\(0\j y \m)\ 



2 _ Q 2 r 2 ^0 -<?™?_, 2^-1 I m even 



2 m - 1 m odd 



Both the relative and absolute intensities are crucially influenced by the parameter c, 
which we have not yet supplied a physical interpretation. To this end, we calculate the 
uncertainty in angular position of the ground state, 

AO 2 = (6 2 - (6) 2 ) 

= c 2 (Q\(a + Jf\0) 
= c 2 (0|(aa f + a f a)|0) 
= c 2 . 

It seems that the parameter which controls the multiple-peak effect is A# = c = y^h/2Iuo, 
that is the extent to which the wavefunction covers the circle. 
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We see that the transitions between the ground and excited states are allowed for all 
excited states in the torsional pendulum, and the distribution of intensities is crucially 
determined by the angular uncertainty In the "stiff pendulum" limit, Uq is large and c< 1. 
The ground state wave function in this case subtends a small angle and motion along the 
periphery of the pendulum is well approximated by the appropriate Cartesian coordinate 
y = rsin.6 ~ r6. The |0) — > |1) transition is by far the strongest, with the other peak 
exponentially suppressed both as a function of c and m. Conversely, in a floppy pendulum, 
the perpendicular motion is important to the response and also the matrix elements to 
higher states become appreciable. The crossover between these limits occurs when the 
angular uncertainty c becomes comparable to 1 radian. 

For values of uncertainty c >1 radian, the maximum intensity is no longer the |0) — > |1) 
transition, but rather the |0) — > |2) begins to dominate the oscillator strength. Figure 121 
shows this in a plot of |(0|j|m)| 2 versus final state quantum number m for several values of 
A6. 

The diagonal ((crf x + af y )/2) optical conductivity is: 

<xiH = ^^\{m\m\ 2 {K^-E m + E Q ) + 5{tiuj + E m -E Q )) 

Vuj I 2 / — ' ml 

m=l 



7r Q r _ c 2 \ m 



e~ c2 —Ac 2 ) 111 1 (5(huj - mhwo) + S(fvjj + mhuj )) 



V 21 m 

m=l 



We can integrate the optical conductivity and find the total oscillator strength: 
Jo V 2/ t^i m - Jo 

oo 

2 m ^ 



71 


Q 2 r 2 


V 


21 


71 


Q 2 r 2 


V 


21 


7T 


Q 2 r 2 


V 


21 


7i 


Q 2 r 2 



Ac 

m 

m=l 



171 ' 2\m-l 



oo 1 

V- 



'c 2 r 



n=0 



= ^1 

V 21 8 
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FIG. 3: The conductivity of the torsional pendulum versus frequency for three values of 
therefore uq). Inset: The same data plotted as a function of final state quantum number. 



We can also find 02(0;) using the Kramer-Kronig relation[4]: 



„(„) = -hp r 



7T 



Uj' — UJ 



m=l ^ 00 



7?) 



! co> 2 — (mwo)" 
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giving for the complex conductivity (er = a\ + 102) 



2 00 / 

X ^— ' m.' 1 



m! a; 2 — (mcun) 

m=l 

Further insight into the physical significance of multiple peaks comes by considering 
the time-dependent current which arises from the application of a short electric field pulse 
E(t) oc 5{t). In this case, E(uS) — E — const, and for t > 0, 

1 f 00 

J(t) = — / e - w V(o;) J B(a;)cL; 



OO 



qAU) p 2 «, m 2 \m-l + 

e y — Ac) cosmt^ot 



2ir 4 ' m 

m=l 

2vr 4 2 l j 

(for i < 0, J(t) = 0, a consequence of causality that is built-in to the Kramers-Kronig 
relations.) This current response is shown for three values of c in Figure |U For small 
values of c, the response of the system is similar to that of a harmonic oscillator, exhibiting 
nearly sinusoidal oscillations for t > 0. Loosening the pendulum (and increasing the angular 
uncertainty A9) effects these dynamics considerably 

While the considerations spelled out here are rather idealized, a reasonable place to seek 
the multiple peak effect in a real system could be the vibrational spectra of nanotubes. When 
a chiral nanotube is formed from a graphene sheet, phonons propagating along the graphene 
lattice basis vectors form a twisting pattern around the nanotube axis. The geometry 
associated with these vibrational degrees of freedom bear similarities to that of Figure ^p. 
It is feasible that another realization of this effect could be found in solids which support 
very soft librational phonon modes. 
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FIG. 4: (a) The current response J(t) to an impulse pulse at time t = 0. (b) The same curves as 
(a), plotted versus the scaled time variable in order to show how the sinusoidal response changes 
as the pendulum is loosened. 



I. APPENDIX 



Here we present some calculational details: 



Jx = ^[H,P X 



^ r rr2 



[L 2 ,cos9] 



2lh 

-(L 2 cos 9 — cos$L 



21 h 

(L 2 cos 9 — L cos 9L + L cos 9L — cos 9U 



2lh 

^j^(L{L, cos 9) + [L, cos9]L). 
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The commutator in this expression is 



d d 
IX, cos #1 = —ih(— cos 9 — cos 9—) 
o9 39 

d d 
= —ih(— sin 9 + cos 9— — cos 9—) 

o9 o9 

= ihsm.9 



and so 



A similar expression follows j y 



^(L[L,cos#] + [L, cos 9]L) 
-^-(Lsm9 + sm9L). 
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